Observation of Two-Neutrino Double-Beta Decay in 136 Xe with EXO-200 
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We report the observation of 

~>21 



two-neutrino double-beta decay in 



°Xe with 



Ti/2=2.11±0.04(stat.)±0.21(sys.) x 10 1 yr. This second order process, predicted by the Standard 
Model, has been observed for several nuclei but not for 136 Xe. The observed decay rate provides 
new input to matrix element calculations and to the search for the more interesting neutrino-less 
double-beta decay, the most sensitive probe for the existence of Majorana particles and the 
measurement of the neutrino mass scale. 



PACS numbers: 23.40.-s, 14.60.Pq 

Several even-even nuclei are stable against ordinary /3 
decay but are unstable for /3/3 decay in which two neu- 
trons are changed into two protons simultaneously. As is 
well known, /3/3 decay can proceed through several modes. 
The allowed process, the two- neutrino mode (2^/3/3), is 
completely described by known physics; its rate was first 
evaluated in Of the other, hypothetical, modes, the 
neutrino-less decay (0^/3/3) is forbidden in the Standard 
Model since it violates conservation of the total lepton 
number. Its observation would constitute proof that 
neutrinos are Majorana leptons Q, unlike all charged 
fermions that are of the Dirac type 0|. Moreover, the 
(V/3/3 decay can proceed only if neutrinos have mass 0] ■ 
Consequently, there is an intense worldwide program of 
experiments aiming at observing the 0v[3/3 mode. The 



relation between the 0^/3/3 half-life and the average Ma- 
jorana neutrino mass requires the evaluation of nuclear 
matrix elements that, while different from those of the 
2i//3/3 decay mode, would benefit from their knowledge. 
Indeed, it has been suggested [H that the theoretical pa- 
rameters contributing to the largest uncertainties in the 
Of/3/3 matrix element calculation can be derived from the 
2^/3/3 decay matrix elements, known once the half-life 
has been experimentally measured. The half-life of the 
2^/3/3 decay depends on details of the nuclear structure 
that are only known approximately [6]. The 2vj3(3 deca; 
has been observed in all important candidate nuclei 
with one notable exception, 136 Xe, which until now had 
only lower limits on the half-life {8j]. The most stringent 
published limit would imply a nuclear matrix element 
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FIG. 1: Drawing of the EXO-200 TPC. The chamber contains 
~175 kg of liquid Xenon enriched to 80.6% in the isotope 136. 



noticeably smaller than those found for other isotopes. 
From an experimental perspective, the 0^/3/3 and 2v/3/3 
modes can be distinguished from the study of the energy 
spectrum of the electrons. The sum energy spectrum is a 
resolution-limited line at the Q- value for Ov/3/3 (2458 keV 
for 136 Xe and a broad continuum for 2is(3f3, so that 
good energy resolution, along with the knowledge of the 
2is/3(3 rate are essential in the search for 0^/3/3. 

The EXO-200 detector, shown in Figure [TJ is a time 
projection chamber [ljj (TPC) using liquid Xe (LXe) 
both as the source of nuclear decays and the detection 
medium. The TPC has the geometry of a cylinder of 
40 cm diameter and 44 cm length, with the drift field 
obtained by biasing a cathode grid dividing the cylinder 
into two identical regions. Each end of the cylinder is 
flared into a conical section, containing two wire grids 
and one array of ~ 250 large-area avalanche photodi- 
odes [ll| (LAAPDs) that allow for simultaneous readout 
of ionization and scintillation in the LXe. Wire grids 
cross at a 60° angle, providing 2-dimensional localization 
and energy readout of each charge deposition. The third 
(longitudinal) coordinate is obtained from the time inter- 
val between the scintillation signal in the LAAPDs and 
the collection of the charge at the grids. A set of field 
shaping rings, lined with reflective teflon tiles, grades the 
field and limits the drift region to two cylinders, each of 
18.3 cm radius and 19.2 cm length. For the data pre- 
sented here the cathode bias was set to -8.0 kV, provid- 
ing a field of 376 V/cm, designed to be uniform to within 
1% over the entire fiducial volume. This low value of 
the electric field provides more stable operation at the 
expense of the ionization energy resolution that is not 



essential for the measurement of the 2^/3/3 mode. 

All components used for the construction of the detec- 
tor were carefully selected for low radioactive content 12 1 
and compatibility with electron drift in LXe. The TPC is 
mounted in the center of a low-background cr yos tat filled 
with - 2400 1 of high-purity HFE7000 fluid [H] serving 
the purpose of innermost radiation shield and heat trans- 
fer fluid. At least 50 cm of HFE7000 (with a density of 
1.8 g/cm 3 at 167 K) separate the TPC from other com- 
ponents. The LXe (and the HFE7000 fluid) is held at 
147 kPa (1100 torr) and 167 K, with possible tempera- 
ture variations < ±0.1 K, by cooling the inner vessel of 
the cryostat with a closed circuit refrigerator. The cryo- 
stat is vacuum insulated and has a total radial thickness 
of 5 cm of low background copper. It is further encased in 
a 25 cm thick low-activity lead shield. Signals from wire 
triplets, spaced 9 mm from each other, and LAAPDs are 
brought out of the cryostat and lead shield, where they 
are amplified, shaped and digitized at 1 MS/s by room 
temperature electronics. The detector infrastructure in- 
cludes a gas phase recirculation system consisting of a 
boiler, a pump, a hot Zr purifier, gas purity monitors 14 1 
and a condenser. A substantial control system maintains 
a very small (< 85 torr) pressure difference across the 
TPC vessel that is built out of 1.37 mm thin copper to 
keep backgrounds low. A calibration system allows the 
insertion of miniaturized radioactive sources to various 
positions immediately outside of the TPC. 

The clean room module housing the cryostat and the 
TPC is surrounded on four sides by an array of 50 mm 
thick plastic scintillator panels [15j. The array detects 
muons traversing the lead shielding with an efficiency 
of 95.9%. EXO-200 is located at a depth of about 
1600 m.w.e. in a salt deposit of the Waste Isolation Pi- 
lot Plant (WIPP), near Carlsbad, NM. The muon flux at 
this site was measured [l6| to be 3.1 x 10~ 7 s _1 cm _2 sr _1 . 
A paper describing the EXO-200 detector in detail is in 
preparation. 

For the data presented here EXO-200 was filled with 
~175 kg of xenon enriched to 80.6 ± 0.1% in the iso- 
tope 136 ( enr Xe). The remaining fraction (19.4%) is the 
isotope 134, the rest of the natural Xe isotopes repre- 
sent negligible contributions. 85 Kr is a radioactive fission 
product with Qp = 687.1 keV and T\/ 2 = 10.8 yr that is 
present in the atmosphere since the nuclear age and gen- 
erally contaminates Xe, as a trace component of natural 
Kr. The EXO-200 enriched Xe was measured 0, El to 
contain (25 ± 3) x 10~ 12 g/g of natural Kr, substantially 
less than the typical concentration of 10~ 8 to 10~ 7 found 
in Xe after distillation from air. 

The data were collected between May 21, 2011 and 
July 9, 2011, for a total of 752.66 hrs of low background 
running. During the same period, about two hours of ev- 
ery day were devoted to detector calibration using 60 Co 
and 228 Th sources. A specific source was inserted each 
day at one out of five "standard" positions near the TPC. 
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FIG. 2: Energy spectrum for a 228 Th calibration source at the 
mid-plane of the TPC, 3 cm outside the LXe volume. The in- 
tensity (vertical scale) is not fit: the agreement between the 
Monte Carlo calculations (solid line) and the data tests the 
accuracy of the simulation against the absolute, National In- 
stitute of Standards and Technology traceable, source activ- 
ity. The energy scale (horizontal axis) has been corrected for 
r e at the time of collection and has had separate energy cal- 
ibrations applied for single- (main plot) and multiple-cluster 
(inset) interactions. 

Typically a source was used to scan these positions over 
a week and then replaced with a different one. Data 
analysis is performed by two independent groups, provid- 
ing cross checks of the results. The detector calibration 
procedure begins by fitting the energy spectra from the 
sources to obtain an electron lifetime (r e ) in the LXe and 
an overall correspondence between the charge and the en- 
ergy deposited in the detector. After an initial phase of 
recirculation of the Xe r e reached 250 /is and it remained 
between 210 and 280 /is in the data set used here (the 
maximum drift time at the field used here is ~ 100 fis). 
T e values are obtained by minimizing the energy resolu- 
tion for source calibration events, occurring at locations 
randomly distributed over the entire LXe volume. The 
dispersion in the r e measurements is incorporated in the 
systematic uncertainty. The daily calibration schedule 
makes it possible to track and correct for changes in r e . 

For this initial analysis only the ionization signal is 
used to measure the energy. The scintillation signals 
recorded by the LAAPDs are used to establish the time 
of the event, identify a particles by their higher light-to- 
charge ratio compared to electron-like events and mea- 
sure a energies. The combined use of scintillation and 
ionization, to obtain the best energy resolution, is under 
development. 

The ability of the TPC to reconstruct energy deposi- 
tions in space is used to remove interactions at the de- 



tector edges where the background is higher. It also dis- 
criminates between single-cluster depositions, character- 
istic of (3(3 and single (3 decays in the bulk of the Xe, from 
multi-cluster ones, generally due to 7-rays that constitute 
the majority of the background. In the present analysis, 
such discrimination only employs one spatial coordinate 
(«15 mm separation) and the time coordinate («17 mm 
separation). The fiducial volume used here contains 63 kg 
of onr Xe (2.26 x 10 26 136 Xe atoms). The detector simula- 
tion, based on GEANT4 , reproduces the energy spec- 
tra taken with calibration sources well. This also applies 
to the single- to multi-cluster assignment obtained with 
the external calibration sources, as illustrated in Figure[2] 
for the case of 228 Th. 

Four full absorption 7 calibration peaks, spanning the 
energy region of interest for this analysis, are derived 
from the 60 Co and 228 Th sources: 1173 keV, 1332 keV, 
2615 keV and 511 keV (annihilation radiation). The 
three high energy 7s provide both single-cluster and 
multiple-cluster event samples. The energy scale is found 
to be slightly (~ 4%) different in the two samples be- 
cause of the non-zero charge collection threshold on in- 
dividual wire triplets. An additional calibration energy 
at 1592 keV is provided by selecting ionization sites pro- 
duced by e + e~ pairs from the highest energy 7s. The 
nature of these energy depositions, however, is different 
from the others, being produced directly by ionization 
in a smaller volume. This type of deposition is analo- 
gous to that expected from (3(3 decay and is found to 
be slightly different from single-cluster depositions from 
7s. This shift is well reproduced by the simulation, once 
induction between neighboring wire triplets and other 
electronics effects are taken into account. After correct- 
ing for these two shifts and the (slowly) time-varying r e 
the energy scale fits well to a linear function. The frac- 
tional residuals from this process are shown in the top 
panel of Figure EH 

The measured energy resolution is <je = 4.5% at 2615 
keV. A parameterization of the resolution function is in- 
corporated into the simulation, as shown in the bottom 
panel of Figure [3J The analysis uses an energy threshold 
of 720 keV, chosen so that both the trigger and event 
reconstruction are fully efficient. Probability Distribu- 
tion Functions (PDFs) for each source and position are 
generated by means of Monte Carlo simulation and com- 
pared to the single- and multiple-cluster data (see Figure 
[2]). This procedure reproduces the activities of the exter- 
nal calibration 7 sources to within ±8% of their known 
activities. 

The data collected during low background running re- 
quires only two selection cuts to remove modest back- 
grounds. Cosmic-ray induced backgrounds are rejected 
by removing events preceded by a veto counter hit in a 
5 ms window. This cut removes 124 events introducing 
a dead time of 0.12%. The decay rate of 222 Rn is in- 
dependently determined to be 4.5±0.5 /iBq kg" 1 from 



4 









• Single-site full energy 

■ Pair production 

a Multi-site full energy 


— 


\ X — — 


e+/e- 




i \ 511 


""Co 


228 Th 








1 








1 








1 




, , , /, , , 1 







350 



500 



1000 1500 2000 2500 
true energy (keV) 



3000 3500 




1000 1500 2000 2500 
true energy (keV) 



3500 




1000 



1500 2000 2500 

reconstructed energy P(3 (keV) 



FIG. 3: Top: fractional residuals between the energy calibra- 
tion points and the linear model discussed in the text. The 
single- (solid line) and multi-cluster (dotted line) uncertainty 
bands are systematic, stemming from the finite accuracy of 
the position reconstruction and the r e correction. The thick 
dashed line represents the central value of the shift predicted 
by the simulation for point-like energy depositions. Bottom: 
measured energy resolution (points) along with a parameter- 
ization (line). 

an a-spectroscopy analysis performed using only scin- 
tillation signals, consistent with f}—a and a-a time co- 
incidence analyses. Similarly, 220 Rn is constrained to 
< 0.04 /zBq kg" 1 (90% CL). In the data set 72 /3-a co- 
incidences are removed. The implementation of this cut 
introduces a 6.3% dead time due to spurious ionization or 
scintillation signals. Events are then classified as single- 
or multi-cluster and energy spectra are obtained for these 
two classes, as shown in FigureHJ The spectra are simul- 
taneously fit to PDFs for the 2vf3(3 decay signal (65% of 
which is above threshold) and various backgrounds using 
an un-binned maximum likelihood method. The 2v{3f3 
PDF is produced using the Fermi function calculation 
given in [2(J. The detector simulation predicts a small 
fraction of the 2vfij3 decay signal to be classified as the 
multi-cluster type because of brehmmstrahlung as well as 
charge collection effects. Background models are devel- 
oped for various components of the detector, inspired by 
screening of materials performed at the time of the de- 
tector's construction and by estimated cosmogenic acti- 
vation. As Figure[4]illustrates, the backgrounds involving 
7 rays are readily identified by their clear multi-cluster 
signature, while the single-cluster spectrum is dominated 
by a large structure with a shape consistent with the 
2vf3ft decay of 136 Xe. The simultaneous likelihood fit to 



FIG. 4: Energy distributions from 752.66 hrs of EXO-200 
single-cluster events (main panel) and multi-cluster events 
(inset). The result of a likelihood fit to a model including 
the 2vf3fi decay and several backgrounds is shown (solid line) 
along with the 2V/3/3 component (shaded region) and some 
prominent background components at the radius of the TPC 
vessel ( 232 Th, long dash; 40 K, dash; 60 Co, dash-dot; 54 Mn, 
thin dash; 65 Zn, thin solid; 238 U chain in equilibrium, dash 
double-dot). Other background components fitting to negligi- 
ble amounts are not shown, for clarity. The energy scale used 
for the main panel is consistent with that of single-cluster, 
/3-like events while the scale of the inset is consistent with the 
multi-cluster events it represents. The combined xV degrees 
of freedom between the model and the data for the two binned 
distributions shown here is 85/90. 



the single and multi-cluster spectra reports a strong sig- 
nal from the 2i//3/3 decay (3886 events) and a dominant 
contamination from 40 K at the location of the TPC ves- 
sel (385 events). Other contributions account for a total 
of less than 650 events, each with a very low significance 
in the fit. These levels of contamination are consistent 
with the material screening measurements (l2| . Taking 
only the single-cluster events into account the single-to- 
background ratio is 9.4 to 1. 

The a-spectroscopy analysis is used to bound any 238 U 
contamination in the bulk LXe. This is important be- 
cause 238 TJ decays are followed (with an average delay of 
~ 35 d) by 234m Pa decays, producing j3s with a Q- value 
of 2195 keV. The a scintillation spectrum is calibrated 
using the lines observed from the 222 Rn chain, obtaining 
a limit for 238 U (and 234m Pa) of < 10 counts for the data 
set shown in Figure HI In addition, a study of the pro- 
duction of fast neutrons resulting in recoils and captures 
in the LXe as well as thermal neutrons resulting in cap- 
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tures is used to bound these backgrounds to < 10 events 
for the data set in the Figure. 

The measured half-life of the 2v(3j3 de- 
cay in 136 Xe obtained by the likelihood fit is 
T 1/2 =2.11±0.04(stat)±0.21(sys)xl0 21 yr, where the 
systematic uncertainty includes contributions from 
the energy calibration (1.8%), multiplicity assignment 
(3.0%), fiducial volume (9.3%) and 7 background models 
(0.6%), added in quadrature. The uncertainty from the 
energy calibration is estimated using a Monte Carlo 
method scanning calibration constants within the range 
illustrated in Figure [3] and re-fitting the spectra, weight- 
ing the fit results by their likelihood value. The same 
method is used to quantify the effect of the multiplicity 
assignment. The fiducial volume uncertainty is deter- 
mined from the fidelity with which calibration events 
are reconstructed within a chosen volume as compared 
to simulation. The 7 background model uncertainty is 
derived from the results of likelihood fits performed with 
a variety of different background hypotheses. 

In Figure E] the fitted values of the 2^/3/3 and the K 
background are shown as functions of the event standoff 
distance from materials other than the LXe (top panel) 
and time in the run. While the 40 K is attenuated by the 
LXe as expected, the 2i//3(3 signal appears to be uniformly 
distributed in the detector and constant in time. 

An exhaustive search for j3 emitters with no 7s, Ti/ 2 > 
2 days and energies of interest yields only two candidates: 
90 Y (supported by 90 Sr) and 188 Re (supported by 188 W). 
It appears a priori unlikely that the bulk of the LXe 
is uniformly contaminated with these isotopes while si- 
multaneously not showing significant evidence for more 
common metallic contaminants such as those from the 
238 U decay chain. Nevertheless, additional test fits are 
performed by incorporating each isotope separately. At 
90% C.L. the 2V/3/3 rate is reduced by less than 7% (30%) 
for the inclusion of 90 Y ( 188 Re). 

In conclusion the initial data taking of EXO-200 has 
provided a clear detection of the 2v/3f3 decay in 136 Xe. 
The measured is significantly lower than the lower 
limits quoted in [8j and translates to a nuclear matrix 
element of 0.019 MeV" 1 , the smallest measured among 
the 2i/f3ft emitters. 
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RFBR in Russia. The collaboration gratefully acknowl- 
edges the hospitality of WIPP. 
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FIG. 5: Top: measured 2u/3j3 decay rate of 136 Xe (large 
points) and largest background contribution ( 40 K, small 
points) as a function of the standoff distance from detector 
components. Bottom: event rates of 2v/3/3 and 40 K decays as 
a function of time. 
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